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SUMMARY 


The purpose of this study was to improve the computational 
efficiency of the impedance- type payload/launch vehicle inter- 
face loads prediction method (Reference 1) developed under tasks 
1 and 5 of Contract NASI -14370 and to evaluate its use as a 
practical design tool. Three goals were defined: 

1) Devise a method of selectively removing from the cal- 
culation process those components of transfer function coeffi- 
cients which do not appreciably affect the interface environment 
in order to make the impedance method operate more efficiently 
in the computer. 

2) Assess the accuracy and convenience of the method for 
determining the effect of design changes (design interactions 
as opposed to completely new designs) in either the payload or 
the support structure. 

3) Investigate the possibility of using the method to 
identify beneficial design changes which may reduce loads on the 
payload. 

In pursuit of these goals, the interface acceleration 
expressions of Equation 28 in Reference 1 were derived by 
direct inversion of the coupled system impedance as expressed 
in terms of booster and payload component modes. From this 
exercise it was concluded that the method will operate more 
efficiently in the computer by expressing the system impedance 
in terms of orthogonal coordinates. Consequently, an approxi- 
mate ortnogonal ization approach was taken whereby eigenvectors 
determined for a baseline design are used to develop uncoupled 
equations of motion, without a new coupled system eigensoluti> 

The formulation is suitable for calculation of dynamic response 
in either the frequency or time domain. Thus, computational 
efficiency was improved while uncertainties involved in component 
mode selection were avoided. 

To assess accuracy, modes of the Shuttle Transportation 
System (STS) landing configuration with two simple payloads 
were used to exercise the method. Modal characteristics of 
the coupled orbiter with baseline payloads were calculated 
by conventional methods. The payload models were perturbated 
and new coupled orbiter and payload modal characteristics were 
calculated by both a conventional method and by using the base- 
line eigensolution to transform from component to approximate 
coupled system mcdes. Structural transfer functions defining 
frequency response of payload loads to unit landing gear 
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sinusoidal input and payload loads time response to unit step 
landing gear input were calculated using both the conventional 
and approximate mode coordinates. Approximate mode results 
compared well with conventional results. 

The approximate modal characteristics were used to calculate 
payload model kinetic and potential energy distributions as de- 
fined in Reference 2. These energy distributions were used to 
identify payload changes which shift system frequencies away 
from those regions identified by transfer function data as bad 
for payload loads. 

It was concluded that: 

1) Analyses integration, manpower, and schedule devoted 

to modal coupling is identical to that required to implement the 
impedance technique equations. Computer costs for the additional 
t step of generating coupled system modes will be recovered by 
decreasing impedance technique cost due to o, prating with diago- 
nal matrices. Even coupled system costs can be avoided for 
model perturbations such as those which occur during design. 

2) The accuracy of the approximate mode impedance technique 
is adequate for determining the effect of design changes. It 

is very convenient to use and works for both statically determi- 
nate and statically indeterminate payload attachments. 

3) Beneficial design changes to reduce loads on payloads 
can be identified by the combined application of impedance tech- 
niques and energy distribution review techniques. 
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INTRODUCTION 


Payload loads analyses are lengthy and expensive tasks. 

Part of the expense can be attributed to intercompany integration 
efforts during a payload design cycle. Customarily, the payload 
contractor generates a dynamic model of his current design and 
transmits that model to the booster contractor. This transmittal 
generally involves at least one working group meeting to identify 
payload nodes, degrees of freedom (dof's), load transformations, 
etc. and to agree on data formats, coupling techniques, computer 
interfacing, etc. This integration effort continues by phone 
and/or mail throughout each load cycle and is repeated for each 
cycle. 

Another part of the expense can be attributed to the number 
of load cycles. Whenever a payload change evolves in the course 
of design activity the effect of the change on payload loads 
usually cannot be assessed without repeating the entire loads 
cycle effort. 

Tasks 1 and 5 of this contract were initiated to develop 
an approach whereby the payload contractor could proceed through 
some portion of his design effort independent of the booster 
contractor and assess the dynamic impact of evolutionary payload 
changes to determine if a repeat of the load cycle is necessary. 

The result of these tasks was the development of an impedance 
expression relating interface accelerations to booster forces 
which can be implemented repeatedly by the payload contractor 
as he sees fit. A method of scaling results from one payload 
model to another on the same booster also was described. Only 
one transmittal of data from the booster contractor to the pay- 
load contractor is required with this approach. This impedance 
technique is reported in Reference 1. 

As in any frequency domain analysis, the impedance technique 
requires the solution of simultaneous linear equations at each 
frequency of interest. For large systems (many dof’s) and broad 
frequency. ranges, these calculations become expensive. One objec- 
tive of the investigation described here was to make the impedance 
method operate more efficiently in the computer. 

An additional objective of this task was to evaluate the 
use of the impedance method as a practical design tool. 

The approach to both objectives is based on the fact that 
computations are less expensive using orthogonal coordinate systems 
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(i. e., diagonal coefficient matrices), and on the assumption 
that a single coupled booster/baseline payload system eigen- 
solution will provide mode shape data adequate for generating 
approximate orthogonal coordinates for all coupled booster/ 
perturbed payloads resulting from evolutionary payload design 
changes. 


4 



SYMBOLS AND ABBREVIATIONS 


dof 

STS 

(hUhMh} 

{W},{F},{h} 

{qMqUq} 

{q},{q},{q} 


degree of freedom 

Shuttle Transportation System 

absolute discrete coordinates 

discrete coordinates relative to a fixed interface 

absolute modal coordinates 

modal coordinates of system constrained at the 
interface 


[4;] matrix of unrestrained component normal modes 

of vibration 

[ 4 >] matrix of restrained component normal modes of 

vibration 

[T] static deflection shapes associated with unit in- 

terface discrete motion 


[I] 

tf] 

[Z] 

M 

W 

W 

[$] 


unity matrix (a diagonal of ones) 

unit matrix with zeros diagonals at interface dof 

Impedance matrix 

Matrix of coupled booster/payload system modes 

diagonal matrix of unrestrained component eigenvalues 

diagonal matrix of restrained component eigenvalues 

diagonal matrix of coupled booster/payload system 
eigenvalues 

a subset of [4;] which couples the booster/revised 
payload component modes 
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V 

[i^] A subset of [vO which does not couple the booster/ 

revised payload component modes. 

PIC Payload Integration Contract 

MMC Martin Marietta Corporation 

[DD] Transformation relating payload loads to absolute 

discrete displacements. 

A defined to be 

^ approximately equal to 

S Laplace operator 

• Therefore 

• • 

Subscripts : 


B booster 

I , interface 

P payload 
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CHOICE OF COORDINATES 


One approach to forming coupled booster/ payload system, 
homogenous, conservative equations of motion starts by de- 
fining the absolute discrete motion, {h}, of all dofs as a 
linear combination of payload/booster interface coordinates, 
{hj}, plus elastic motion relative to that interface,^}. 

That is: 

{h} = [T] (hj } +{ R } (1) 

where [T] defines the static deflection shapes associated with 
unit interface displacements (e. g. a Gi yan reduction trans- 
formation). Equation (1) may be partitioned into booster, B, 



Equation (2) is then used, along with the respective mass, 
[M], and stiffness, [K], matrices to form the equations of motion 
by classical energy methods. The results appear as, 




At this point, several options exist for the choice of normal 
mode generation. Frequency, the following choice is made: 

1) Constrain {h.} and {h-l and solve the resulting eigenproblem 
described by the (3,3) partition of Equation (3) to obtain eigen- 
values (££], and eigenvectors, [? p J, of the payload constrained 

at the interface. 
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2) Delete T p M p T p and T p K p T p . Then, constrain fh p } and sol''e 
the resulting eigenprbblem described by the (1,1), (1,2), (2,1), 
and (2,2) partitions of Equation (3) to obtain eigenvalues, 
[io^],and eigenvectors, [<j>g], of the unrestrained booster without 

payload effects. 

3) From the results of steps 1 and 2 form 


h 


^ IB 

— 


q B 

V 

h 

r 

r 

I 



q P 


4) With proper normalization, use Equation (4) to transform 
Equation (3) into 


Ih ^b t p m p t p^ib 

^pVp ~P 
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(5) 


In practice, step 1 above is performed by the payload contractor 
who transmits his results to the booster contractor each time a 
payload design is perturbed. The booster contractor then performs 
steps 3 and 4. If booster design does not change, some savings 
in intercompany communication and data transmittal can be made by 
transmitting booster data (step 2) from booster to payload contractor 
or.j time only. However, the effort required to accomplish steps 
3 :nd 4 is unchanged. 


Now component damping and booster forces can be added to Equation 
/ 5) and Laplace transformation into the frequency domain (with 
».aro initial conditions) yields, 


S“I+S2-,g +< P IB Tp(S" , Mp+Kp)Tpd) IB 

S 2 <PpIpM p T p 4) 1B 


S ^P 

S z I+S2<'a* 2 +(o p 


q B (s) 

q p (s) 


(F(S)l 


( 6 ) 


where S is the Laplace operate. 



A solution for the modal accelerations can be obtained by 
direct inversion of the acceleration coefficient matrix. Inver 
sion by parts (Reference 4) allows the result. 


{q B (S)> = S 2 [(S 2 I+S2? B w B +^)^J B Tj(S 2 M p +Kp)T p « IB 


-1 


"1 


- s \'{b t p m p^p <? p (s2i+S2c pV U) p ) ^pi p Mp T p 4»i B D *[«J]{F B (S)} 


By defining, 

[ IMP(S ) 1 - (LTp] T (S Z [M p HK p ])[TpHTjM p Ip$ p ] 

(S 2 fS 2 I+S2^ p +w p 3 WjlpMpTp] 

and manipulating Equation (7), the form 

[I-S 2 (S 2 I+S2c B co B +o3 B ) 1 <j>{ B [IMP(S)]<J, IB ]{q B (S)) 

- S 2 [S 2 I+S2;o)^u) 2 ] [4>J]{F b (S)} 

can be obtained. If Equation (9) is premultiplied by [4 >t B ] 
and definitions 

(hj(S)} = [<J>jg]{q B (S)} 

[PADM(S)] = [4. IB ](S 2 fS 2 I+S2i; B u) B +^ )[<f>{ B ] 

[TADM(S)] = [^ IB ](S 2 CS 2 I+S2c B to B +co B g )[<f>J E ] 

are made. Equation (9) becomes identical to the impedance 
expression given by Equation (28) of Reference 1 which was 
developed in Task 1 and 5. That is. 



(13) 


[I-PADH(S)IMP(S)]{h I (S)} = [TADM(S)]{F(S)> 

Thus, the effort required to develop the coupled equations of 
motion is no more than that required to set up the impedance 
technique. The form of Equation (9), and its equivalent in 
Reference 1, is due to the choice of coordinates. That is, 
component normal mode coordinates were used, rather than coupled 
system normal mode coordinates, in order to save the costs of 
the coupled system eigenproblem, as defined by Equation (5). 

If the eigenproblem of Equation (5) is solved and the result- 
ing eigenvectors ^re denoted by corresponding eigenvalues 
are denoted by and coupled system normal mode coordinates 
are denoted by {£}. Transformation of Equation (5) to system 
mode coordinates with proper normalization of [♦], and adding 
damping and booster forces and performing the Laplace Trans- 
formation, as was done to form Equation (6), yields, 



f(S 2 I+S2cfl+fi 2 )3{C(S)} = |> T ] 



(F b (S)} 


The expression for discrete interface accelerations is obtained 
from Equation (14) by inverting a diagonal matrix. That is. 


(14) 


fhj(S)} 


S 2 [4> BI 0][i^]f(S 2 +S2?n+fi 2 ) J|> T ] 



(F B (S)} 


The comparable solution in component normal mode coordinates is 
obtained by multiplying Equation (7) by [<t> B j]* 


(15) 


{hjCS)} = S 2 C4> BI 3[(S 2 I+S2c B ^ B +ojg)+(bj B Tj(S 2 M p +K p )T p 4> I8 

“5 4* jgl"pMp!p ( { | p(S I + S2CpWp+Wp) ^pIpMpTpifjg] 


uJ](f b (s)> 


(16) 


The advantages of Equation (15) over Equation (16) with respect 
to computation cost is obvious. Cost avoidance associated with 
not solving the coupled sysetm eigenproblem is rapidly over- 
shadowed by repeated inversions of a coupled impedance matrix. 
Because of this, the original task approach for improvement of 
computational efficiency was to identify component modes which 
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do not contribute significantly to the results and delete them, 
thereby allowing cheaper inversions of a smaller matrix. The 
approach was abandoned for the following reasons: 

1) Frequency truncation of component modes above some 
cutoff value (e.g., 50 Hz), as is customarily done, usually does 
not decrease the problem to an economical size. 


2) Although some booster component modes within the range 
of interest can be deleted due to their negligible amplitude at 

the point of excitation, (i.e., [<j>J RFgl-O in the n^ mode) 


i r 


significant size reduction by this approach is not likely. 


3) Since payload changes are causing payload component 
frequency changes throughout the design period a payload component 
mode judged insignificant in one design could become a dominant 
contribution if a change caused a frequency shift which greatly 
coupled that mode to a retained booster component mode. 

Consequently, the approach was revised to that of identifying 
an economical approximate eigensolution suitable for estimating 
gross effects of evolutionary design changes. 
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APPROXIMATE COORDINATES 


If one assumes that Equation (5), in the previous section, 
was developed for the baseline payload model at the beginning of 
a design effort, and that [i|>] and are the eigenvectors and 
eigenvalues respectively calculated for that baseline system by 
conventional methods, then transformation of the equations of 
motion into these coupled system normal mode coordinates gives 


[i ](€>+ 


= U ,T ] 


Fb> 


(17) 


The computational advantages of this diagonal system have been 
discussed in the previous section. 


If, in the course of design, the payload changes somehow, 
the payload model will change and new matrices to describe payload 
mass, stiffness, constrained component modes and frequencies, and 
static deflection shapes associated with interface motion will be 
formed. If these new matrices are denoted by subscript 2, the 
coupled system equations appear as 


I+ ^IB T P 2 M P 2 T P 2 <{> IB 

(J ’ib t p 2 m p 2 i p 2 

i’p 

V Z 

T ^ 

d>p Ip Mp Tp $td 

L 2 V Z V Z *2 1 


i 



“eWb^^Jp^ib 

0 

_2 

k i 


^1 

0 

(jL)n 

4 

n 1 

1 

0 

- 


He 

b p. 


08 ) 


A conventional eigensolution from the homogeneous form of 
Equation (18) can be calculated and an imnediate transformation to 
the diagonal form of Equation (17) made. (This will be done to 
obtain comparative results). 


Or, if the same number of booster and payload component 
modes are used to define the revised coupled system as were used 
in the baseline coupled system the substitution of 



(19) 
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into Equation (18) can be made. If all modes of the baseline 
system were calculated and included in Equation (19), the result- 
ing equations of motion 



i+4, ib t f 2 m p 2 t p 2 ^ib 


IT “ 

4, IB T P 2 M P 2 I P 2 ^P 2 





Mp Tp ip 
2 v 2 v 2 v 2 


IB 



+ 



2, ,T t T „ 

W B 4> IB T P 2 K P 2 ^IB 



(20) 


would have no loss of generality. Adding a damping matrix 



w j 

to Equation (20) and defining the coupled mass and stiffness 

ma + rices as[ H r ] and[ K r Respectively, Equation (20) can be 
l 2 l 2 

written in the form 



Equation (21) can be used to develop a set of approximately 
orthogonal coordinates to use in lieu of solving the revised 
system eigenproblem. This is possible because, although [t|] 
are not the true eigenvectors of the perturbed system, they are 
not fa** from it. If [<[] are "in the neighborhood" of the eigen- 
vect s of the perturbed system, 

»!“*■[ % 

(Reference 3) given by the ratios of the diagonal terms in the 
stiffness matrix divided by the corresponding diagonal terms in 


a good estimate of the eigen- 
can be obtained from the Rayleigh quotients 


the mass matrix as indicated in Equation (21^. If the coupling 
produced by [’Jin Equation (21) is sufficiently small ,a good 
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estimate of that system of equations is given by 


-1 


-1 


[!]{£> + fMeqJ [2c c Meqn c }{£} + [MeqJ [Me^ }{£;} 


* CMeof [«' B T ][4. B T ]{F B } (22) 

where [Meq] represents the diagonals of M r ^], [MeqQ r J represents 

T l 2 l 2 
the diagonals of [y K r ^]and f2£ r Meqn r ^represents the diagonals 

of[* T C r *]. 


There is no criteria for determining what is sufficiently 
small or whether or not the vectors of the baseline system are in 
the neighborhood of the eigenvectors of the revised system. But, 
when one considers that successful engineering decisions are often 
based on imperfect orthogonality (in both analysis and test) one 
might expect that use of the baseline system coupled modes as 
if they were orthogonal with respect to the revised system mass 
and stiffness will produce useful intermediate payload design 
information. The remainder of this report generates an example 
problem to illustrate that this expectation can be satisfied 
even when payload revisions are fairly large. 
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DYNAMIC MODEL 


Many theories have been demonstrated by using very simple 
dynamic models. Frequently, the examples are trivial, failing to 
provide a valid test. To avoid such a trap, an existing landing 
configuration STS set of modal characteristics will be used as 
the booster model and two payloads will be indeterminately con- 
nected to the cargo bay. Although the payloads are fictitious 
and lacking great detail, they interface with the STS in the same 
way as true payloads and are representative of their mass and 
frequency content. 

All STS lander model data was obtained from Martin Marietta 
Corporation (MMC) Denver Division in-house computer files which 
were established under the Payload Integration Contract (PIC) for 
the United States Air Force. The data was handled in a manner 
similar to that used on PIC. That is, all modes in the frequency 
range 20 Hz and below were included as normal modes. All modes 
above 20 Hz were included as 14 residual flexibility modes 
(Reference 5) associated with 14 interface dof to which the two 
payloads attach. On the PIC Project, modes up to and including 
50 Hz were included as normal modes and the remainder as residual 
flexibility modes. Budget restrictions here would not allow use 
of the large size coupled system that results from this frequency 
cutoff, but the 20 Hz cutoff that was used serves well to illustrate 
the point. Table A-l lists 84 STS lander model frequencies, the 
last 14 of which represent residual flexibility modes. Table A-2 
identifies cargo bay and landing gear degrees of freedom in the 
STS lander model. 

Figures 1 and 2 describe the geometry and physical properties 
used to generate the two baseline payload models. These models 
were attached to the orbiter by inertial coupling with elastic 
properties described in terms of constrained (at the STS inter- 
face) modes. Equations (2) and (4) define this coupling. Stand- 
ard MMC methodology and computer codes were used to generate the 
payload models and to calculate their modal characteristics. 

Tables B-l and B-2 define the payload dof tables. Tables B-3 and B-4 
record the modal characteristics of the baseline payload models. 
Tables B-5 and B-C define the constraint modes ([T] in the equations) 
for the two payloads as obtained from the respective payload 
stiffness matrices by static reduction techniques. 

The STS lander and baseline payload models were coupled by 
the method previously described and modal characteristics of the 
coupled system were calculated. Table 1 lists the resulting 
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FIGURE 2: Description of Baseline 

Payload Model B 
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frequencies under the column headed "COUPLED MODES." Table 1 also 
identifies the component modes which contribute more than 5% of 
the kinetic energy of each coupled mode, under columns headed 
"MAJOR CONTRIBUTING COMPONENT MODES." Subheadings under this 
identify contributing component modes as "IDENTIFICATION," % of 
contribution as "CONTRIB," and uncoupled component frequency as 
"FREQ." STS lander normal modes are identified as "ORBXX," where 
XX = 1,50. STS lander residual flexibility modes are identified 
as "ITERFYY" where YY = 1,14. Payload constrained modes are 
identified as "AMODZZ" or "BMODZZ" where ZZ = 1,10. Thus, the 
total number of component modes and coupled system eigenproblem 
size was 84. The result of interest from the baseline coupled 
system are the eigenvectors X which transform component 
mode coordinates to coupled system mode coordinates. The 
matrix is often referred to as "modal modes." This 
84 x 84 matrix is presented as Table C-l. The product f 
[Mq][v \ where[ Mj;]denotes the coupled baseline system 
results in 1.0 on all diagonals with the largest off dia 
element being 10' 13 . The product [ij> ] t [Kc][^ 1 ,where[ K £.j jci.otes 
the coupled baseline system stiffness, results in eigenvalues on 
all diagonals with the largest off diagonal element being 10 -3 . 

These results will be compared later to orthogonality checks of 
[ 4 j ]with respect to the coupled perturbed system mass and stiff- 
ness matrices. Note for now the relationship between the energy 
contribution in Table 1 and the relative magnitude of terms in 
tv’] as shown in Table C-l. Where a coupled mode is shown in 
Table 1 to consist primarily of a single component, the column 
in[y] corresponding to that mode has a value near 1. in the 
row corresponding to the component mode. For example, element 
7,12 of v is 0.9941 while the 12th coupled mode is shown to be 
99 . 093 % orbiter mode 7 (0RB7) which is the 7th component mode. 

Even when coupling does occur, as in couoled mode 15, only those 
rows of [ip] corresponding to the participating component modes 
contain significant numbers. 

To obtain the perturbed payload models, arbitrary changes 
of mass and stiffness were made. Figures 3 and 4 describe the 
perturbed models. Degree of freedom tables are unchanged from 
the baseline descriptions given by Tables B-l and B-2. Modal charac- 
teristics of the perturbed payloads are presented in Tables D-l 
and D-2,and constraint modes in Tables D-3 and D-4. Note that 
perturbed model constraint modes differ from the baseline due 
to stiffness changes in the indetermina trt load path to the 
interface. 

The STS lander and perturbed payload models were coupled by 
the same method used for the baseline payload. This coupled 
perturbed system was handled in two ways. 
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FIGURE 3: Description of Pertjrbed 

Payload Model A 
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FIGURE 4: Description of Perturbed 

Payload Model B 



1) Modal characteristics were calculated by the same method 
as that used on the baseline system. Resulting frequencies are 
listed in Table 2, and "modal modes" are recorded as Table E-l. 
These results provide the true orthogonal coordinate system for 
the response analyses comparison which will be discussed in the 
dynamic model excitation section. Many of the coupled modes are 
seen to be primarily booster modes. This result is consistent 
with our PIC results where approximately half the coupled system 
modes are predominantly booster modes. 


2} Orthogonality checks, [s^][M r ][^] and [**][ K r ][i] 
where [M r ]and [K r ]are the coupled perturbed system mass and 

2 4 

stiffness matrices respectively, were formed. They are presented 
as Tables F-l and F-2. Rayleigh Quotients were formed from the 


diagonals of these matrices (i.e., 
were converted to frequency (i.e.. 


Kp /M r , 

2jj 4jj 



j = 1,84). They 




) 


and correlated by mode shape with the true eigensolution results 
from step 1. The correlation is presented in Table 3, where the 
Rayleigh Quotients are superimposed on the table of Perturbed 
Coupled System Frequencies. This good correlation verifies that 
the coupled system "modal modes" obtained from the baseline pay- 
loads eigensolution provides a very good estimate of the coupled 
system "modal modes" for the perturbed system. 
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DYNAMIC MODEL EXCITATION 


With feasibility of the aporoximate coordinates approach 
established by the preceeding review of the perturbed dynamic 
model, both time and frequency responses were calculated using 
first the true orthogonal coordinates (step 1 above) and then 
repeated with the proposed approximate coordinates (step 2 above). 

Interface loads were selected as the variable to be calcula- 
ted and landing gear forces were selected as the excitation source 
Since only comparative results were desired, a unit force was 
applied at all frequencies in frequency response calculations. 
Forces were applied in the vertical direction at all three struts. 
No attempt was made to represent any realistic landing gear forces 
In general, the equations of motion were set up in the form, 

CI](U+[2^]{C}4^J{U=[^ T ]M {F bl } (23) 

and the loads equation were set up as 

{L} = [DD][T<J, BI ilpMU) (24) 

where [DD] is a transformation relating payload loads to absolute 
discrete displacement. 

Loads in Element 3 of both perturbed payloads (see Figures 3 and 
4) were selected for presentation because their major axes coin- 
cide with the model axes, making it easier to correlate load 
magnitudes with the excitation forces. Figures 5 through 16 
present frequency responses of these loads. Figures 17 through 
23 present time responses of these loads. The "a" subscript 
on these figures denotes that the true orthogonal coordinates 
were used. The "b" subscript denotes that the approximate coordi- 
nate were used. Table 4 decodes the six load names used in 
Figures 5 through 28. 

Since input forces were arbitrary only relative magnitude of 
loads can be considered when reviewing results. In general, loads 
obta.ned using orthogonal coordinates agree fairly well with 
corresponding loads obtained using approximate coordinates, parti- 
cularly in the direction of excitation where loads are the 
largest. Further improvement of frequency response results may 
be possible by using a smaller frequency increment in the analysis 
(.127 HZ was used in the interest of economy). 



Table 4. Notation for Response Data Figures 


LOAD NAME CODES 

PX = Axial Load in X direction 
MTX = Moment about X axis 
VY = Shear in y direction 
MTZ = Moment about Z axis 
VZ = Shear in Z direction 
MTY = Moment about Y axis 

(See Figures 1 and 2 for description 
of axes) 

COORDINATE SYSTEM DESIGNATIONS 
"a" denotes true orthogonal coordinates 
"b" denotes approximate coordinates 
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Figure 5b. P/L A, PX, Frequency Response 








loads leoeeo ortho srs p/l a i tv loads freq resp 
Figure 61. P/L A, MTX, Frequency Response 
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Figure 7a. P/L A, VY, Frequency response 
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Figure 7b. P/L A, V Y, Frequency Respon 
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Fioure 8a. P/L A, MTZ, Frequency Response 
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F inure 8b. P/L A, MTZ, Frequency Response 










Figure 9a. P/L A, VZ. Frequency Response 






Figure 9b. P/L A, VZ, Frequency Response 






Fimirp 10a. P/L A. f'TY, Freauency Response 
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Figure 10b. P/L A, f!TY, Frequency Response 
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Figure 11a. P/L B, PX, Frequency Response 
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Figure 11b. P/L B, PX, Frequency Response 




47 


Figure 12a. P/L B, HTX, Frequency Response 







Figure 12b. P/L B, KTX, Frequency Response 
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Figure 13a. P/L B, VY, Frequency Respot se 





Ficjure 13b. P/L B, VY, Frequency Response 
















Figure 14a. P/L B, MTZ, Frequency Response 
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Figure 14b. P/L B.MTZ, Frequency Response 
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Figure 15a. P/L B, VX, Fiequency Response 





Figure 15b. P/L B, VX, Frequency Response 
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Figure 16a. P/L B, MTY, Frequency Response 
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Figure 16b. P/L B, I'.TY, Frequency Response 
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Figure 18a. P/L A, f-TTX, Tine Response 


700*00 



60 


LOADS 160E0O APROX SYS P/L A I /F LOADS RFSP TO STEP 
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Figure 19b. P/L A, VY, Time Response 
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Figure 20a. P/L A, MTZ.Tirne Response 
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Figure 22a. P/L A, liTY, Time Response 
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LOADS I6OE0O APBOX SYS P/ L 0 I /r LOADS RESP TO STEP 
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Figure 24a. P/L B, HTX, Time Response 
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Figure 25a. P/L B, VY, Time Response 
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LOADS 16OE80 ORTHO SYS P/ L B |/F LOADS RESP TO STEP 
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Figure 28a. P/L B, MTY, Time Response 
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IDENTIFICATION OF BENEFICIAL MODEL CHANGES 


If the frequency content of the excitation source is known 
and transfer functions from a critical load to the point of exci- 
tation is calculated, it is easy to determine what frequency 
changes in the model would be beneficial to load reduction. 

The problem then, is to identify what model changes will pro- 
duce the beneficial frequency shift. A solution to this problem 
has been suggested in Reference 2. The approach is conveniently 
applicaDle during payload design activity. 


Frequency shifts are produced by mass and/or stiffness changes. 
In aerospace dp;igns, minimum mass is almost always a goal for 
performance rposes. Therefore, adding mass to produce a fre- 
quency shift probably would not be allowed. Removing mass probably 
would not be po>sible. Therefore, further discussion will 
concentrate on the identification of beneficial stiffness changes. 

(A similar discussion of mass changes is possible, however.) 


The homogeneous form of equation (17) provides a convenient 
place to start this discussion. [ft 2 .] in this case was formed Dy 




I B T P K P T P^ I B 


— 2 


+ [oJ]&P*p] (25) 

The explanation is easier if a statically determinate interface 
is assumed for the moment. With this assumption, the j th diago- 
nal term of Equation (25) can be written as 


o 2 

'“'j 




l^pJppKVj 


(26) 


If the jth couDled mode is primarily a booster mode, the right 
hand side of Equation (26) will be small. Where payload parti- 
cipation in the jth coupled mode is important the right hand 
side of Equation (26) approaches ft?. In some cases 
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( 27 ) 



In this case, if examination of excitation source and payload 
transfer functions indicate a shift of the jth coupled mode is 
beneficial, a shift in the ith payload mode frequency is obvious- 
ly in order. However, w 2 p; also affects the frequency of other 
coupled system modes. An estimate of the overall affect of chancing 
payload frequencies can be obtained as follows. 

A A 

Define [ft 2 ] and [ y] as the new coupled system modal charac- 
teristics obtained as the result of perturbing each of the n 
payload components frequencies by 6.or-, where i = l,n. Then 
the new form of Equation (26) is 1 P 


(28) 


As was shown in the discussion of the Dynamic Model, the unperturbed 
coupled system eigenvectors provide a good estimate for the per- 
turbed system. Therefore, Equation (28) can be written as. 


[4][5pSp[*p]^a I ]-[4][u»|][» B )-[4]C|]['!.p3't4n 2 ] 


(29) 


where [Aft 2 ] are the difference between the desired coupled eigen- 
values and those given by the current model. Since the diagonal 
terms of the left hand side of Equation (29) can be written as 


L p i j ^ 


: i J 1 6 Pi u Pi’ g ' 


2 

Pi j 


i-vpiJjt 


wp](5 p }, 


(30) 


an approximate solution for the component mode frequency scaling 
can be obtained from 


on 

[y 2 ] is of size j x i where is the number of coupled modes and 

i is the number of payload component modes and j>i. A subset of 
rows can be selected from [*„..] to obtain a square nonsingular 

matrix, allowing a solution of { 6p} . A least squares solution, 

involving all data in [bp.] also is possible. This possibility 

is discussed in Reverence 2. 
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A similar approach to determine the payload modifications 
required to achieve the approximate payload frequency hanges 
also is described in detail in Reference 2. 


The matrix [typ.] is presented as Table 5. It was formed 
using the baseline^doupled system modes, ([4<] has been 

presented previously as Table C-l). The circled row numbers 
in Table 5 were selected to form the non-singular matrix re- 
quired in Equation 31. The corresponding [wt] matrix of Equa- 
tion 31 was assembled from the baseline payload data. These 
baseline model frequencies are summarized as column 2 of 
Table 6. The {Air} required for Equation 31 were calculated 
from the difference between the baseline and perturbed coupled 
system eigenvalues. Column 5 of Table 6 presents the baseline 
payload mode scaling factors, { 6 }, as predicted by Equation 31 
to produce the perturbed coupled system frequencies. Column 4 
of Table 6 presents the actual baseline payload frequency per- 
turbation factors that resulted from changing from the baseline 
to perturbeJ component models. The small differences between 
the predicted factor (column 5) and the actual factor (column 4) 
shows that the approach quite accurately identifies component 
model changes required to produce a desired coupled system 


frequency shift. 
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Table 5 . Coefficients for payload component frequency scaling 
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Table 6. Payload Component Frequency 
Scaling Comparison 
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2) obtained by solving Equation (31) for {5 



CONCLUSIONS 


1. The choice or Q n orthogonal coordinate system will allow the 
impedance method to operate more efficiently in the computer. 

2. The eigenproblem required to determine ai. orthogonal coordinate 
system can be formed as a by-product of the impedance method 
set up effort. 

3. If one solution for the coupled boostar/payload system is 
obtained, the resulting "modal modes" can be used to form 
approximate orthogonal coordinate systems very economically 
for each payload perturbation that occurs during the design 
process. 

4. The "modal modes" can be used to determine payload component 
mode frequency changes that are beneficial to payload loads. 

5. The method works for statically indeterminate payload attach- 
ments as well as determinate ones. 
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APPENDIX B 


BASELINE PAYLOAD MODEL DATA 
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10.33. IB CLOCK TIME 

MOOL3 RUN WITH LANOER ANO TWO FAKE PAYLOADS 37.431 SIC. CPTIME 

USE FOR IMPEDANCE STUDY MODELA AND MODE LB 30993 SEC. PPTIMC 
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Table B-5. Baseline payload model A constraint mode 
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APPENDIX C 


BASELINE COUPLED SYSTEM MODAL MODES 
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11.33.32 CLOCK TIME 

MOOL3 RUN WITH LANOER AND TWO FAKE PAYLOADS 149.610 SEC. CPTIME 

USE FOR IMPEDANCE STUDY MOOELA ANO MODELS 48783 SEC. PPTIME 
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11-33.33 CLOCK TIME 

M00L3 RUN WITH LANOER ANO TWO FAKE PAYLOADS 149.971 SEC. CPTIME 

USE FOR IMPEDANCE STUOV MODELA AND MOOELB 48831 SEC PPTIME 
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APPENDIX D 


PERTURBED PAYLOAD MODEL DATA 



Table D-l . Perturbed payload model A modal characteristics 
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19.48.49 CLOCK TIME 

WOOL 3 RUN WITH LANDER AND TWO FAKE PAYLOADS 161.326 SEC. CPTJME 

USE FOR IWPEDANCE STUDY MODE L A6 AND MODE LB6 40937 SEC. PPTIME 
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19.40.49 CLOCK TIME 

HOD 13 HUN WITH LANDER AND TWO FAKE RAY LOADS 161.071 SEC. CRT I ME 

USE FOR IMPEDANCE STUDY MODE LA6 AND MODE LB6 40909 SEC PR T I ME 
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APPENDIX F 

CROSS ORTHOGONALITY CHECKS 



Table F-l . Check orthogonality of baseline modal modes 
with respect to perturbed system mass matrix 
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NEW FR'IO AND MOOES FRO* SELECTED MOOES 28.243 SEC. CPTIME 

FORCE COEFFICIENTS FORMED 5669 SEC. PPTIME 


8 

? 



t 

o o o o o 5 

nt©©®®©© 

oooooooo 

W Appp«a 

OOOOOOOO 

^^©©©©©o 

OOOOOOOO 

«©»^©©©TOTO 

oooopoco 

® r> © © 

o o o o 

w 

U WWWIUIU 

triwOOF 
O CHO t I"* ^ 

o * m r> tr 

yj u u uj m m m ^ 
•'*»TOt-«©TO® 

*-nhQ(7i>»'T 

O — r-00©®n 

Mi Mi Ui Ui Ui M| m| mi 
© ct — ®ctTOftf- 

-er>Qe>©m«« 

«®OnOt®« 

w w u uj uj uj uj uj 
'*’» r'nO0»©<T' 
tnt'-o*- to 
r- n © © r» 

lU UJ UI UJ !U uj uj uj 
TOOOQTOO nr-© 
ncPOvsns 
h0NPO0O*' 

© ui Ui ui 
V O 0 N 
— — cm r- 

TO — — Oi 


*- ' c n s o 

nnc P** wnoi 

w-wntwnp 

fWN«t®N« 

«®0j — — r*© — 

- - TO © 


lO (0 TO © ® © 

OOOOOO 

oooooooo 

tototo©to©to© 

OOO O O O O C 

»ntifi»pps 

OOOOOOOO 

P®PAhPPl8 

oooooooo 

n © © © 
O O OO 

01 

« * o — o © 

0 ** t - ® o 

Mi m mi mi m ^ Mi Mi 

©"* Dno®p- 
— TO TO © TO 0» TO — 

p^^On-PO 

rtnO® no® O 

TO ® TO to n — — © 
TOOt^nctoTOft 

Mi ui Ui Ui uj uj uj uj 
ppr-ct(TPjit 
©r- D«« © 0 n 
n©^©©(7»©'® 

i ii iii 1. 1 m uj uj uj uj 
OOO i — Q©TOO 
O— OTOOrt — r- 
® — nr>r-TOO® 

Ui © © © 

O ® ® 

© — TO TO 
A fi A B 


* *"f>0no 

vnntDnnh 

•-« — ©to©— r» 

— — © — ©o© 

©ex — nncxn® 

- O' © — 


t«# tn* 

OOOOOO 

^nppp^np 

OOOOOOOO 

©®©©©toto© 

OOOOOOOO 

vr*^©©^©© 

OOOOOOOO 

® ® ® © TO © ® r* 

oooooooo 

® ® © © 

o O D O 

• 

uuitu WUMD 
ro p- r> r> cm <ji 
. o — — r- r> — 
© m f“ <0 ct 

— ©©—©•—« 
pm-PO'pn® 

®<*©£>®m — o 

Mi mi mi iij Mi m m Mi 

0*O©totoo«o 

t^rOn^fiN 
ct © © — n to r« © 

©<*« — ©« — © 
©no-® — ®©to 

UJ UJ Ui LlI UJ IU UJ U i 
O — TO©CtOi© — 
0O^sn0nn 
— ©r'©o»©C® 

© © © © 
— O « TO 
© n n 

TO O 0» <T> 


^ M 

r'etc-©TOTO-*“« 

©r-0»TOTOTO — — 

— — ro — o» o» i . « 

— — — nnr-©n 

— — no 


« OC000 

OOOOOO 

OOOOOOOO 

oooooooo 

®©w©©tr©© 

OOOOOOOO 

©TO®©©©©© 

oooooooo 

® ® ® TO 

o o o O 

r* 

m m |y |^j m m 

lO0f>t OO 
® n « « <0 n 

~ ® m t-* ® 

Mi Mi Mi Ui Mi Ui U> Mi 

POoOnnp- 

PPMTPNnNP 

Mi Mi uJ Ui Ul Mi UI 

r~ m r> O' ® — © © 

©OTOTOTOTOTOr- 

PMhClpn®»- 

Mi Mi Ui (M UJ Ui UJ Ui 

r-nr* — dj^toO 
to©-®©©- ©n 
r'®<^(T»oi — oi — 

kJUJJ. UJUiUlLiJuJ 

0» — ©«)nQ®c® 
POAnPOPc* 
®TO0®0®’^‘P» 

UJ U © Ui 

TO r- O O 
® O ® n 
r* to © r- 




n — — ct«ctr~ct 

*-®®ci®nn- 

— ®?-TO 

TO - — *- 


* i0 to i0 © © 

OOOOOO 

® * ** p ^ r. p © 

OOoOOOOO 

to©totototoc-to 

OOOOOOOO 

©©©©©©©© 

OOOOOOOO 

©©©TO©©TOr~ 

OOOOOOOO 

© ® © TO 

O O O O 

to 

uj uj u m uj m 
0O^n«t 

00D«OC 

yj u Ui UJ UJ U UJ UJ 
**t*'r"ipor>®i0 

■,t <n - n - o f' 

t r C p p « 

Mi |y Mi UJ Ml Mi Mi Uj 

p®p«® - no 
n®nnp® n t 
O0>©O«h-r'O 

m UJ Ui UJ UJ UJ UJ UJ 
r* ® (7> © co ® ®r- 
0»AOn*-w0 
r^c*®®ot©CM« 

ui uj ul Ul uJ Ui UJ u 1 
0>0>TOl s ' — TOOtru 
© cv r- ® CD ® 

— cxr-nn© — — 

© © © © 
© n to r- 

<N © W 0> 
none® 


— ft <0 10 IT - 

— Ct — O'- — 

m;*) — — ^^«TO 


— TO— — — OiCtTO 

e — ® *- 


ip to »- 10 r* % 
OOOOOC 

t P» » « P P C) «5 
OOOOOOOO 

tototo®©©toto 

O O OjO O O O o 

®C»©®©©’CI^ 

OOOOOOCO 

©®©©TOQ®TO 

OOOOOOOO 

® ® © © 

o o c o 

© 

UJUMUUIUIU 

ct © to c 0> w 

< ij. v <ji ^ 

cr o» to ® O 

UJ UJ Ui UJ UJ kJ UJ u> 
BNpr. p®o>(* 
noipn^t&pn 
rttfOffl^PO® 

yj UJ uj uj uj uj UJ V 

©©©cMQnnuj 
— .T©®© — noo 

«oo®o®^^ 

Uj Ui Ui Mi Mi UJ UJ UJ 
f'~00®® — r-tP 
Onrtft«®n 
0nntPP«O 

u ui ui Ui uj uj Ui uj 
c*o«®TO-“Qr‘0% 

— ©r-P'np®® 
r^w®O^0®C» 

Ui UJ © © 

to — © r« 
— n n — 
© Cm O r- 


^ rt *- p « 

® o ft t- nno*- 

n ® * n n n p n 

« — ©©© — ©I s - 

© — — O'D — Oto 

— — — a) 


ip tr 0 *i © n © 
O O o O O O O 

ooooooooo 

tBnttPlDPB 

OOOOOOOOO 

©©n®-®-©©® 

OOOOOOOOO 

®-®®©©©r'©TO 
OOOOOOO^ ^ 

® n ® ® 

O O O O 

® 

w®t OP»** 
r- TO n o ® ft n 
n o a> o r- ® tP 

ui uj uj w uj uj uj Ui w 
n®«p®0n®O 
rtnOwPOOW® 
fit POPPn®'' 

U W UJ UJ UJ Ui IMJ UJ UI 

®®0O — ©OCDK 
®nNC v®apn 

i s -©c*®<j'cr> — ©to 
© rt p* © cm ff- n - © 
n ft od s p T'ct-toc! 

Uj uj w uj uj o u 
r- to o> cm r- - 0» 

0O0A®©pN .. 
O*~(TJ0'©©©*‘ n 

© © © UJ 

® O r* (n 
n ® O to 
© p> C n 



— — © — n — © ct o 

— a — — — ©.-ct — 

«©n® — ©ct— m 

-0P)*PJ'-0Mt 

® - n - 

Q 

VPhpPAf 

OCOOOOC 

ppnnppppp 

OCOOOOOOO 

to©toto©Q©©© 

OOOOOOOOO 

©©^®©©©©ao 

OOOOOOOOO 

©®n® TO © TO © 
OOOOOOOOO 

© ©TO © 

o O O O 

is 

2 

Ul KM Mi Mi UJ yj Mj 

vvnnhO^ 

MCDPO>*- flf 

— ^ * © G ® © 

Mi Mi Ui ui Ml/ Mi yj Mi 

Cv«pncthp* 
pciono® »P*- 
n*’-rO’"riPe 

Mi Mi Ui UJ Mi UJ ui Ui UJ 

r» — to— r-Q©On 
TOTOO"~0>Qooct© 
P»«ppon so 

UJ UJ UJ UJ UJ Ui Ui Ui UJ 

®«© — ©r-©r'Q 

oo — n©ni-r<o 
'“©©n©©n©© 

Ml Ui Mi Ml Ui UJ VM Mi UJ 

«B0O0nOn*- 
Oto — r-r>— to — <T> 
p 0 n n* t O ^ t- 

W Ul W UJ 
TO O CD n 

© n rt 
TO D TO 

O 

o 

t CtPPW^ 0 

o»c?ft — m©TO® — 

— ct® — — — nr-- 

i ■ 

0jr't»©©©ctO0' 

® — — — CMOJ — — O' 

® ® (C co 


m to © ® w ip © 

ooooooo 

ttntP’-ttn 

OOOOOOOOO 

©TOTOn©©®©® 

OOOOOOOOO 

©©^»©r^®©©® 

OOOOOOOOO 

p n t © to © © to © 
OOOOOOOOO 

® ® © © 

o o o o 

— ■ 
O ft 

a 

►- w 

O9«n««o 
.1 <n r» n r- <* 
*■ O ® 10 ct to 

vnnofippnn 
O O' » ff> © t pmi 
0iO*“O^C»conrj 

i^j i. I ill yj ||J II i 1. 1 yj y 

' r* © — O to O' — ® 
®©T“n«®TO® 
© © O' - — ^- © Oj 

r'-AjQfc'ji — ®©©n 

n© — 0©C®cr> — 
© n © cn ® — 

ui mi m Mi yj Mi UJ Ui UJ 

C't^t^r^r-TOfT®tT> 
®0(P®r-TO®nO 
r^ — o*ct— ©r-®® 

© UJ © UJ 

O TO ® TO 

— 1 - TO n 
TO TO TO - 

o 

N 

• tit 


© — ^*“*©*-r- * 

• t« i « i 

r- ct a* « ®r — © © © 

OJ — — » CMMnoi® 

ft C' r- ® 


V 0PQf>C0 

O O O O 0 o o 

nnt ®p t npp 
OOOOOOOOO 

v^v©©tnv©tf> 

OOOCOOOCO 

®®®-©©©®©© 

OOOOOOOOO 

©©®TOTO©®TOr> 

ooooooooo 

© to n to 

o o o o 

® - 
m 

111 Ui U 141 1*1 H) Ul 

® *t ® C o O «- 

r* ® to o ip «*- to 
o n r- o to © © 

Mi Ml Mi 881 Mi Ui UJ u* Ui 
«np0ppppp 
nnOPOti - u ® 
n w p « t — © © rt 

UlVkUiJJU WUUJUI 

— ® — r'nooc>f-'0t 
®O^OnPft««- 

©®©O^O<T>0‘^ 
Qr-n — non®® 
O f TO ® n TO © n 

M; UJ UJ Ui Mi Ui 1 i UJ UJ 

no; — oo®0®TOp 
— ®n — r-®o&©0 
ici/ttttONMr 

© © © © 
© O' © n 
ft n co r- 
O ft TO N- 

z 

n to ®* *• to « to 

i« t i 

*- ft to a> to *- c r> ~ 

®nrirttpf*«® 

®© — ©nn© — © 

iN©ct — r. — ©ctr- 

TO to — ft 

- 

n ro to © to r» to 

►orntpch# 

-ftn^©«r^« 

«n® p»h ® 

**wr>®©TO r-« 

^ ft n 

o 

o 

© © © © r © 

nrtnwcicinp'w 
© © © © ^ p P P P 

wnr >r)nr>nn 
©©©©©©©©© 

©©©©©©©©© 

©©©©»©©©© 

©©©©©©©TO© 

TO TO © TO 
© ©TO © 


17. 54,23 CLOCK TIME 
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17.54.23 CLOCK TIME 

HV* F»EO A NO MOOES FDOM SELECTED MODES 2* . 993 SEC . CPTIME 

FORCE COEFFICIENTS FOflMIO 6769 SEC. PPTIME 
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17.54.23 CLOCK TIME 

FREO AND WOOES FROM SELECTED MODES 29.354 SEC. CPTIME 

FORCE COEFFICIENTS FORMED 5817 SEC. PPTIME 
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17.54 24 CLOCK TIME 

NEW FREQ A NO MOOES FROM SELECTED MOOES 29.687 SEC. CPTIME 

FORCE COEFFICIENTS FORMED 5817 SEC. PPTIME 
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le F-2. Check orthogonality of baseline modal modes 

with respect to perturbed system stiffness matrix 
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17 . 54.24 CLOCK TIME 

FREO ANO MOOES FROM SELECTED MOOES 30.739 SEC. CPTIME 

FORCE COEFFICIENTS FORMED 5965 SEC. PPTIME 
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BUN NO. OBBB3B Table F-2. (Continued) PAaE N0 - 28 

17.64.J7 CLOCK TINE 

NEK FRfO 4 NO MOOES FROM SELECTED MOOES 34.414 SEC. CPTIME 

FORCE COEFFICIENTS FORMED 6373 SEC. PPTIME 
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0RBR3R Table F-2. (Continued) PAOC NO. 20 

*7.84.27 CLOCK TIME 

NIK FREQ AND MOOES FROM SELECTED MODES 34.782 SEC. CPTIME 

FORCE COEFFICIENTS FORMED 842* SEC. PRTIME 
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